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Abstract 

We study theoretically the detection and possible utilization of electric current-induced mechan- 
ical torques in ferromagnetic-normal metal heterostructures generated by spin-flip scattering or the 
absorption of transverse spin currents by a ferromagnet. To this end, we analyze the DC voltage 
signals over a spin valve driven by an AC current. In agreement with recent studies, this "recti- 
fication", measured as a function of AC frequency and applied magnetic field, contains important 
information on the magnetostatics and -dynamics. Subsequently, we show that the vibrations ex- 
cited by spin-transfer to the lattice can be detected as a splitting of the DC voltage resonance. 
Finally, we propose a concept for a spin-transfer-driven electric nanomotor based on integrating 
metallic nanowires with carbon nanotubes, in which the current-induced torques generate a rotary 
motion 

PACS numbers: 76.50.+g, 72.25.Ba, 85.85.+j 
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I. INTRODUCTION 



An electromotor is an apparatus that generates rotational motion with electrical cur- 
rents, as demonstrated by Faraday in 1821. Under miniaturization, the torques generated 
by 0rsted magnetic fields scale unfavorably compared to the increased friction, rendering 
the Faraday motor an unfeasible concept at the nanoscale. Alternative concept for nanoma- 
chines are based on electrostatic forces,^ thermal fiuctuations,^, torques induced by circularly 
polarized light'^ and angular momentum transfer by spin-polarized currents." In this article, 
we elaborate on the last idea. 

A spin-polarized current carries an angular momentum current: 



where P = (/^ — /^) //is the polarization of the charge current I = + 1^, e is the electron 
charge and Tq = Ih/(2e). If transferred completely to the lattice, mechanical torques are 
created of the same magnitude, which can be relatively large in nanoscale structures. Since 
spin currents are routinely excited in magnetoelectronic devices such as spin valves, we pose 
here the question whether current-induced mechanical torques can be detected and utilized 
in such structures. We conclude that the resonant magnetomechanical coupling studied 
earlier-i^ should indeed be observable in spin valve structures, paving the way for applications 
such as high frequency actuators and transducers of mechanical motion. Furthermore, we 
propose a design for a spin-transfer driven electric nanomotor based on carbon nanotubes.- 

The resonant rectification of a current in spin valve as a function of an applied AC 
frequency has been found experimentally to form a rich source of information about the 
magnetization dynamics in spin valve structures.-i^ Kupferschmidt et al. found theoretically 
that the spin-pumping by the magnetization dynamics^^ significantly modify these spectra.— 
We suggest that the AC-DC conversion in spin valves can be used to detect vibrational modes 
excited by the spin-polarized currents. 

In this manuscript, we address the theory of spin valves excited by AC currents, show 
how to include the effects of the magnetovibrational coupling, and predict signatures of the 
current-induced mechanical torques. We also share our ideas how these torques could drive a 
rotary (rather than vibrational) motion, i.e. an electric nanomotor. The manuscript is orga- 
nized as follows. In Sec. II, we calculate (position-dependent) mechanical torques generated 
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by the spin-flip dissipation of a spin current injected by a ferromagnet into a normal metal. 
In Sec. Ill, we study the spin-transfer mechanical torques resulting from the absorption of 
transverse spin-currents by a ferromagnet. We suggest to employ the recently reported diode 
effect^i^ in F(erromagnet)|N(ormal)|F(erromagnet) metal spin valves to detect the vibra- 
tions created by spin-transfer mechanical torques. We calculate the nonlinear DC voltage 
emanating from the spin-transfer driven ferromagnetic resonance with and without the res- 
onant magnetovibrational coupling. Spin transfer torques non-collinear with planes formed 
by principal axis of anisotropics can deform the resonant line shape of DC voltage as a 
function of the AC current bias frequency. In Appendix B and C, we calculate the extra 
DC voltage caused by spin-pumping and prove that it can play a key role in sufficiently thin 
films. In Appendix B, we also conclude that the magnetovibrational coupling is observable 
by virtue of the spin-pumping even in asymmetric N|F|N heterostructures. In Sec. IV, we 
propose a spin-transfer driven nanomotor concept based on integrating metallic nanowires 
with carbon nanotubes. 



II. MECHANICAL TORQUES DUE TO DISSIPATION OF SPIN CURRENTS 

Consider a normal-metal diffusive wire or nanostructured pillar into which a spin accu- 
mulation has been injected via an electrically biased ferromagnetic contact (Fig. [I]). Spin 
Is and charge Jq currents can be conveniently related via 2 x 2-matrices in Pauli spin space 
/ = (l/o + a ■ Is)/2, where a is the vector or Pauli spin matrices and 1 is the 2x2 unit 
matrix. Spin-orbit interactions or magnetic impurities cause spin-flip scattering that can 
be parametrized by a spin-flip relaxation time tsf- the local (vector) spin accumulation 
is related to the spin current density = Is/{eS) (e is the electron charge and 5* is the 
cross-section of the wire) by the angular momentum conservation law: 

dt^-'^ dyAf~ tsf- ^ ^ 

The dissipated angular momentum per unit length, r = {h/2) SAfUg/tgf, where A/" is the 
density of states at the Fermi-level, is transferred as a mechanical torque to the lattice. 
In the conflguration sketched in Fig. [H the injected spin accumulation \fis\ = fJ's and the 
mechanical torque are polarized in the y-direction. Newton's Law for the mechanical motion 
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of the substrate then reads^^ 

5^ = + ^s^^^^ (3) 

ot-' oy^ 2 tsf 

where ^{y, t) is the angle of torsion, / = Ix+Iz {iz = f x^dzdx, Iz = f z^dzdx) is the moment 
of inertia of the cross-section at y relative to its center of mass, p is the mass density, C is 
an elastic constant defined by the shape and material of the wire (C = for a circular 

cross section with radius R, /i is the Lame constant). Eq. ([3]) must be complemented by the 
boundary conditions = 0) at the clamping points. 

We now concentrate on a bimetal wire consisting of a ferromagnetic and a normal metal 
(Fig. [T]). We assume for simplicity here that the bulk resistances of the wires are much larger 



than the interface resistance (Ref. llJ shows how interfaces can be taken into account) to the 
extent that we may disregard the latter. The magnetization and mechanical motion is much 
slower than the relaxation scattering, therefore, we can consider only the parametrically 
stationary limit. The charge current density jo is conserved {dyjo = 0) and Eq. ((2]) reduces 
to 

dy-^' tsf hS- ^ ^ 

In the normal metal, charge and spin currents are governed by Pick's Laws = Af^ D^dyfi^ 
and = D^dynf respectively, where D'^ is the diffusion constant and the index 
indicates the normal metal, leading to the diffusion equations: 

— A/'^D^;/^-n —\fNr,N..N_KrNN,.N_^ 
Qy2-'^ ^ /^O -U, Q^^J^ fJ's --^V p^lt^f- 

where rjy is the mechanical torque per unit length for the normal metal. In a ferro- 
magnet (F), the particle and spin currents are = {N'^ D^dyfi^ + D^dyP^ /2 and 
jf = Yadyi^M^ jj,^ — Aff fii)/2, where -D^^) is the diffusion constants for spin-up (- 
down) electrons and J^^^-^ is the corresponding spin-up (-down) density of states. The 
diffusion equation in a ferromagnet then reads: 

leading to 

— L>^(/it -/ij = (/it-/ij/tf^ = r^/ 2*^-^^ 

where Tp is the mechanical torques per unit length for the ferromagnetic metal, 
2D^DlMp/ {AffD^ + AffDf) and Af^ = (Aff + Aff)/2. 
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Figure 1: Spin transfer from spin-polarized currents to the lattice in a lieterostructure consisting 
of two (ferromagnetic and normal metal) sections connected to mechanically clamped normal and 
ferromagnetic metal reservoirs. Polarized currents lead to torques twisting the middle section. 

We can now calculate the mechanical torques by solving the diffusion equations, requiring 
the continuity of distribution functions and conservation of the spin and charge currents at 
the interface.—!^ We also require no spin accumulation at the connection with reservoirs. 
For the ferromagnetic {y < 0) and normal {y > 0) metals we find, respectively: 

Psinh [{LF + y)/Q/smh{LF/0 



/f, coth(L^//fJ + vC coth(L^//fJ ' 
i;Psinh [{Ln - y)/C] /smh{Lr,/li 



tn (y) = Tf 



(5) 



(6) 



'''If, coth(L^//f,) + < coth(L^//fJ 
where P = [G^ — Gi) / + G^) is the current polarization of the ferromagnet here defined 
in terms of the spin-up and spin-down conductances G^ and G^, if. 



DHf^ and C = 

D^t^j are the spin-diffusion lengths in the ferromagnet and normal metals, respectively, 
and Tq has been introduced in Eq. Lp and L^v are the lengths of the ferromagnet and 
normal metals respectively, and v = [-^Ntfj] / [-^Ftff] ■ For the configuration sketched in 
Fig. [H the mechanical torque is directed along the y axis. The torque density is discontinuous 
when V ^ 1. As shown in Fig. 121 the mechanical torques are enhanced close to the F|N 
interface on the scale defined by the spin-diffusion length in both ferromagnet and normal 
metals. The mechanical torques change sign with the electric currents direction. 
Integrating Eqs. (15|6|) leads to the total torque acting on the wire: 



Tp = Tn 



P/frf tanh(L^/2/ 



sd) 



vPlZtanh{LM/2a 



/f,coth(L^/Q+<coth(L^// 



sd) 



(7) 
(8) 




Figure 2: The mechanical torque per unit length along the F|N wire; P = 0.7, = lOnm, 
/f^ = 20nm, Lp/l'^d = 3, Ln/I^^ = 5, and v = 2. 

When Lp ^ if^ and ^ we obtain Tsp = Tp + Tj^ = Ptq as expected by the complete 
dissipation of the spin current in this limit (Eq. ([T])). By ultrasensitive displacement detec- 
tion, it should be possible to observe the mechanical strain caused by the spin-flip torques^ 
in the setup of Fig. [TJ 

III. GENERATION AND DETECTION OF THE MECHANICAL TORQUES DUE 
TO SPIN TRANSFER 

In the previous section, we studied mechanical torques arising from spin-flip relaxation 
processes within the bulk of the metals. In contrast, in this section, we will consider struc- 
tures (see Fig. |3]) in which the spin transfer is dominated by dephasing processes (leading to 
absorption of transversely polarized currents)- at the normal metal [ferromagnetic interfaces, 
whereas the spin-flip processes in the bulk materials are disregarded (which can be repaired 
easily if necessary). Such a device is superior to the wires of the previous section in generat- 
ing mechanical torques when the structures are smaller than the spin-diffusion length. We 
propose here to induce and detect the magneto-vibrational modes^ driven by spin-transfer 
torques in devices such as shown in Fig. |3l ie. a doubly clamped heterostructure in the 
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Figure 3: Model sample to study spin transfer to the lattice consisting of ferromagnets (Fl and 
F2) and NM, i.e.. tunnel junctions or normal metal). Fl and the normal metal N are electron 
reservoirs that are mechanically clamped. The vector diagram shows the reference frames used in 
calculations with respect to the magnetizations. The magnetization of Fl is fixed, Mi = M^^^-g^, 
and that of F2 is variable, M2 = M. 



shape of a bar with a ferromagnetic load in the center. In our set-up, the mechanical torque 
is generated in the ferromagnet F2 by the transversely polarized spin current (ferromagnet 
Fl is supposed to be clamped by the substrate). The absorbed spin-angular momentum is 
transferred to the lattice by the magnetic shape and crystal anisotropics. In the regime of 
the resonant magnetovibrational coupling, the mechanical motion is strongly affected by the 
magnetization dynamics created by spin-transfer torques,- which, in turn affect the trans- 
port properties. We predict that such vibrations should be detectable by current-driven 
ferromagnetic resonance experiments reported in Refs. 9 and 10. We start by analyzing the 
DC voltage signals over a spin valve resulting from the "rectification" of an AC current by 
the precessing magnetization (without magnetovibrational coupling, e.g. setups from Refs. 



y 



and 



lOl ) as an intermediate step. A general analysis of the DC signals in the regime of the 



resonant magnetovibrational coupling is subsequently given. We also estimate the maximum 
torques that can be created in these structures. 
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The derivation below is carried out for a system in which the ferromagnet F2 is inserted 
into the normal metal island. However, the rectification of the AC current and the excitation 
of magnetovibrational modes also occur (and can be treated by our methods) in a system in 
which the ferromagnet is attached to the side or on top of the normal metal wire ("spin-fiip 
transistor") 



A. Ferromagnetic resonance driven by spin-transfer torques 

Initially, we assume that the mechanical modes are not excited in the device depicted in 
Fig. |3l One (the"hard") layer magnetization of Fl is assumed completely fixed by shape or 
crystal anisotropics to the direction Mi = Mfixcd with unit vector mfixcd- We are interested 
in the dynamics of the magnetization direction m = M2/M5 of the middle ("soft") layer F2, 
where IM2I = Mg is the constant saturation magnetization, in the presence of an AC electric 
current bias. We consider the soft layer to be at resonance,^- but note that the soft and hard 
layers can exchange their roles as a function of frequency {e.g. in Fig. |5]). We model the 
magnetization dynamics by the generalized Landau-Lifshitz-Gilbert (LLC) equation in the 
macrospin approximation that is augmented by the current-dependent spin-transfer torque:— 

dM a fdM\ h lit) 



dt 



, . a , . f dM.\ h lit) , , , , ^^ . ^ 

-7MxHeff+— Mx ( — 1+7^-^ [^iHi X (mfixed X m) +7^2 (infixed X m)] , (9) 



where lit) is the time dependent current through the system, r/i describes the efficiency of 
conventional spin transfer, and 772 parametrizes an "effective spin-transfer exchange field". 
When the spacer is an insulator, as in the experiment of the Tsukuba group,- the parameter 
r/i is a constant governed by the expressions such as derived by Slonczewski.— The "effective 
field" 772 has been less well investigated in magnetic tunnel junctions and is treated here as 
an adjustable parameter. When the spacer is a normal metal in a configuration sketched in 
Fig. |3l the LLC equation for the soft layer reads: 



dt 



.,M X H., + X [^'j - X [(I., + 1„) x m] (10) 



where the spin currents leaving the soft ferromagnet can be computed by magnetoelectronic 
circuit theory^^ as 

/i(2) = (Gt + GJ(/if ) - ^^l^'^) + {G, - G,)(Mf ) - /x^(^)) • m (11) 
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- (2m X ^^(2)) X j^^^ _ x /x.^^^)) ^.^ (12) 

where a is the Gilbert damping constant, 7 is the gyromagnetic ratio, and G^ describe 
the conventional spin-dependent conductances limited by bulk and interface scattering and 
G^l = Gr + iGi is the interface mixing conductance of the ferromagnet, /io and are the 
chemical potential and spin-accumulation in the normal metals, respectively, the spin Isi(2) 
and charge /i(2) currents correspond to the current (spin) flow into the normal metal node 
1(2). It is shown in Appendix A that Eqs. and (II Op are equivalent when we allow for an 
angle dependence of the parameters ^71(2)- The efficiencies of the spin-transfer torque rji{6)^ 
and the "effective spin-transfer field" ^72(6') mainly depend on the real and imaginary part of 
the mixing conductance, respectively. 

Choosing the z'-axis along the equilibrium direction (Mq) of the soft layer and the z'-axis 
perpendicular to both magnetizations (see Fig. we expand the free energy close to the 
equilibrium direction of the magnetization as 

F(M) = F(Mo) + N^,Ml,/2 + Ny,Ml,/2 + N^,y,M,,My, (13) 

such that the effective magnetic field 

Heff = -dF/dM = - {N,,M,. + N,,y.My.)^ - {Ny.My. + N,,y,M,,)y' (14) 

where N^i, Ny/, and N^iyi are parameters characterizing the energy of the macrospin ex- 
citations. The free energy Eq. f lT3l) is not diagonal in the basis of the magnetizations 
along the x' and y' axis, but can be diagonalized in another basis rotated by the angle 
(f) around the z'-axis, so that N^' = N^cos'^(f) + NySin^(f), Ny' = Ny cos'^ (j) + N^sm^ (j), 
Nx'y' = {N^ — A^^) cos sin 0. Here, N^ and A^^ are the eigenvalues of the expansion tensor 
in Eq. ffT^ . We consider here an arbitrary direction of an external magnetic field Hq. In 
general, the components of Hq contribute to Eq. (1131) in a non-trivial way. As an illus- 
tration let us consider an external field along the z-axis, which can be included as follows: 
N^, = A^o + Ho /Ms and Ny> = N^, + Hq/Ms, where A^^^^^^ are elements of the magnetic 
anisotropy tensor. 

The linear response of the magnetization to an AC current perturbation is given by the 
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response functions Xx'i = (Mxi/I)^ and Xy'i = {My'/I)i^. From Eqs. (pfT 

% sin 61 ir]2UJ + jMsT^, . . 

_ sin 9 -ir]iuj + -fMsTy> 



■'m 



where Kn is the volume of the magnet, = -^'^N^N^M'^ = -f'^M^{N^>Ny> - N^,y,), = 
r]i{Ny> + N^'y'a) + r]2{N^'y' - Ny>a), Ty> = r]i{N^,y, + A^^;'") + miN^' - N^'y'Ct), and the 
damping parameter is modified by the anisotropics as 

, {N^ + N^)/2 {N,, + Ny,)/2 

a = a — = a (17) 

^JN^Ny' ^jN.'Ny. - Nl^, 

From Eqs. (I15lll6p . we conclude that both the field effect 772, characteristic for tunnel 
junctions,^ and non-collinear anisotropics Nx'y' cause a phase shift in the magnetization 
response Xy/i(a;) that is relevant to the rectification effect, as demonstrated below. At res- 
onance, the parameters Y^' and Yyi are the components of the out-of-phase magnetization 
response. 

The resulting magnetization dynamics causes oscillations of the magnetoresistance 
i?(mi(t), m2(t)) = R{cos6{t)) of the multilayer structure (details of the calculations are 
given in Appendix A). In the presence of an AC current bias I{t) = Iq Re e*'^*, the resistance 
has an oscillating component that in the linear approximation reads: 

i?(cos e{t)) ^ R{iy) - sin A e{t) (18) 

= R{u) - sme^^myit) = R{u) - A gin^^^M Re(e-*Xy'i), (19) 
leading to a nonlinear phase-sensitive effect in the voltage across the sample 

U = R{t)I{t) = R{u)I{t) - ^^|M(Joe-*Xy'i + loe-'-'X'Me'-' + /oe"'-*) (20) 

^P,Re[{l + e^'-')Xy'i], (21) 

where the parameter u = cos 6 describes the equilibrium configuration of the magnetizations 
and 6{t) = 6 + A9(t) (to lowest order AO ^ rriyi). The magnetization dynamics is thus 
manifest in the nonlinear response, i.e. the zero and second harmonic components of the 
voltage across the sample: 



_ lysine dRju) T2-2n^^i^^ 7 

2Ms Of ou 2e2MsVmUJm 



{uj - Wm)2 + a'2u;^ [uj - Lo^y + a''^uj_ 



(22) 
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lysine dR{iy) u^^u^^ 2 ■ 2 n^^i^) ^ 7 

= —— — -— ixy'iMI = /qsiii e- 



As pointed out in Ref. 0, the DC voltage Uq can be interpreted as a diode-like rectifica- 
tion. The amplitudes Uq and U2uj show a resonant enhancement close to Wm (note that f/g 
corresponds to Vmix in Ref. ilO and its sign corresponds to a current flow from soft to hard 
ferromagnetic layer). Eq. (122|) is a linear combination of the symmetric and antisymmetric 
Lorentzians (see Fig. Hj). By fitting the DC voltage Uq{u)) to a linear combination of the two 
curves in Fig. HI we can determine the parameter Ty/ ^ 'riiNx>y' + r]2Nx' that is affected by the 
effective field 772 and the non-coUinear anisotropy N^'y'. In Fig. [5l we demonstrate how the 
resonance becomes skewed merely by the non-collinear shape anisotropy (the spin-transfer 
torque does not lie in a plane formed by principal axes of the shape anisotropy). The res- 
onant layer with the magnetization M2 corresponds to the harder layer in this plot (as it 
was discussed in Sec. IIIA, the harder and softer layers can switch their roles with a proper 
choice of the AC current frequency). The opposite scenario,- without anisotropics, the effec- 
tive field r]2 can still cause an antisymmetric Lorentzian signal since Ty/ ^ ri2Nx' = ri2HQ/Ms, 



as reported by Ref. 



35 



B. Magnetovibrational resonance driven and detected by spin-transfer torques 

In the following, we repeat the derivation of the previous section but allow for a torsional 
degree of freedom along the current flow. This requires a free-standing conducting structure, 
such as a doubly clamped cantilever. We already discussed the coupled dynamics of the 
magnetization and the lattice for a small magnet at the free tip of a cantilever that is 
clamped on the other side.^ However, here we consider the limit in which the magnet is 
heavy compared to the cantilever. This reduces the mechanical resonance frequencies, but 
increases the magnetomechanical coupling strength. The normal metal lattice serves as a 
spring and the ferromagnet is the load. We first consider a well-aligned structure for which 
the primed and unprimed coordinate systems in Fig. E] coincide. A mechanical torsion 
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(co-co )/(a'co ) 



Figure 4: The dimensionless DC voltage over the sample biased by an AC current has two compo- 



nents that are plotted as a function of AC current frequency; Vq = Iq sin 
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diy 2e 2M,K, a'uj„ 



profile described by the angle (p{y) increases the free energy as: 

F(M) = F(Mo) + Kn (iVx [M. + MM^)? /2 + NyMy/2 

+N,y[M, + MMO)]My) + ^ j (^^^ dy, (24) 

where C is an elastic constant defined by the shape of the cross-section and the material 
of the normal metal link (C = for a long plate with thickness a much smaller than 

width d, a d, fi is the Lame constant for the normal metal), v^(0) is the torsion angle 
at the middle magnetic section. The coefficients A^^^, A^^ and A^^^j, have the same meaning 
as in the previous section (one expects Nr^y = in setup of Fig. [6l however, we keep N^y 
in order to use our results for more general configurations). The width L' of the magnetic 
layer along the axis y is supposed to be small compared to the length of the normal metal 
links L(^ L'), so internal strains and deformations in the magnetic section are disregarded. 
The integration is therefore carried out from one clamping point y = —L/2 to the other at 
y = L/2, excluding the ferromagnetic layer. With Eq. (l24l) . the Landau-Lifshitz-Gilbert 
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Figure 5: The resonance peaks are skewed when the spin-transfer torque does not He in a plane 
formed by principal axes of the shape anisotropy. We illustrate this by exciting the magnetization 
M2 of the harder (resonant) layer with an easy plane anisotropy by AC currents, while the soft 
layer magnetization is forced into the directions = : 7r/20 : vr/lO : 7r/4 . by an external magnetic 
field. Here r/2 = 0. 



equation has to be modified as: 
dM , ^ a , , fdM 

-7M X Heff + — M X 



dt 



dt 



+ 7 



h i{t) 



[rjiin X (infixed x m) + ?72(mfixed x m) 



Heff = - (AT.M, + N^yMy + Ar,M,¥p(0)) x 

- {NyMy + N^yM, + N,yM,^{Q)) J, 



(25) 

(26) 
(27) 



where the derivative {dNl/dt)^ is defined in the reference frame of the magnet, since most 
proposed mechanisms for Gilbert damping act on the magnetization motion relative to the 
underlying lattice only. This can be taken into account by 

'dM\ dM rf(^(0) 



dt 



dt dt 



(28) 



The equation of mechanical torsional motion of the normal metal strip is> 

dy^ ^ dt^ ' 
13 



.12 



(29) 




Figure 6: A device to detect magnetovibrational coupling caused by magnetic form and crystal 
anisotropics. An external magnetic field Hq is used to tune the FMR frequency. 

where / = J(2;^ + x'^)dzdx ~ ad^ /12 is again the moment of inertia of the cross-section 
about its center of mass and p is the mass density. The oscillating solution has the form = 
{Aism.{ky) + cos(/i;?/))e*'^*, where k = uj/c is the wave number, c = 2cx,a/d = \JC /{pi), 
and Ct = \/ 1^1 p is the sound velocity of the transverse mode. The free constants and A2 
are determined by the boundary conditions. 

The first condition v^|y=L/2(-L/2) = corresponds to perfect clamping at the boundaries. 
The second one, at the connection with the ferromagnetic load, can be obtained from the 
variational principle applied to Eq. (p^ . and corresponds to the torque Cdip/dy\y=±Li 12 
exerted by the magnetization at the interface to the normal metal link: 

where Ip is the moment of inertia of the magnetic load (for a thin plate of mass M, we can 
approximate it as L'Ippp = Md'j^/ 12, where dp is its width, see Fig. E]), and d!f/dy\y=±L'/2 
is the derivative of the torsion angle of the normal metal link at the connection with the 
ferromagnet. This boundary condition is equivalent to the conservation of the mechanical 
angular momentum written for the magnetic load. With Cd(f/dy\y=±L'/2 = Ck(f cot kL, we 
obtain: 

Ck^iO) cot(fcL) - L'lppp'^^ = 1 + X Ho^ |„ (31) 

where the right hand side describes the magnetically induced torques that to leading order 
equal the Gilbert and anisotropy torques. In the absence of the magnetic anisotropics and 
Gilbert damping, the two terms dNl/dt and 7M x Hq cancel each other. 

By only considering the left hand side of Eq. ( 13T|) . we can find the mechanical resonance 
frequency as = \/ Ck cot{kL) / (L'Ippp) ^ ^JC /{LL'Ippp), where the approximation 
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holds when uj^ is smaller than the resonant frequency of the normal metal link with a free 
end (which corresponds to the heavy-load limit). In frequency space, the left hand side 
of Eq. f l3ip can be expressed in terms of the dimensionless response function F{uj) of the 
mechanical subsystem to an oscillating torque of frequency u applied to the load: 

F{u) = ul/ [{u'' - ul + 2iPu)\ (32) 

where /3 is a phenomenological damping constant describing dissipation in Eq. ( 129|) and it 
is related to the quality factor Q of oscillator at the resonance frequency as Q = Ue/ (2/3) 
(at 1 GHz Q ~ 500).— In terms of F{uj), the mechanical response function reads: 

{ip/T)^ = {l/Ck cot{kL)) F{u) ^ {L/C)F{u), (33) 

where T is the torque externally exerted on the load. 

To first order in the mechanical and magnetization oscillations, Eqs. (125|) and (13T|) lead 
us to the following response function: 



2eVra 

-i7]iuj + r^7 [Ms + g (Ho/N,) Fjuj)] 
lJ^ — iJ^ + ^ia'ijjijj^ + kFiuj) 



(34) 
(35) 



where 



and 



2eVm — + 2ia'ujujTn + AF{uj) ' 
A = g[u^-HouJiN,Ms)], (36) 

g = M^VmN, {l/Ckcot{kL)) (37) 
« M,V„iV,L/C = N,{L/ay{VjVi){M^/fi), 

Vm and Vi are the volumes of the magnetic load and the normal metal spring, respectively. 
L and a are the largest and the smallest dimensions of the normal metal links, k = Ue/c. We 



recover here the results from Ref. 



• 2, i 



implying that the spin polarized current is equivalent 



to an external rf field along the x-axis applied to a magnetic cantilever. Since the ratio 
Vm/Vi can be made much larger compared to the limit of light load considered in Ref. O, 
we conclude that by making the normal metal links thinner (in setup of Ref. 6 we have to 
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make the cantilever thinner) we reduce the stiffness of the device, which results in a better 
sensitivity and stronger coupling. The reduced stiffness leads to a drop in the resonance 
frequency that in turn can be compensated by making the structures smaller. 

The nonlinear response to an AC current can now be found by substituting Eq. 
into the first parts of Eqs. fl22ll23p . We conclude that the most conspicuous feature of the 
magnetovibrational coupling is the formation of a magnetopolariton and the splitting of the 
ferromagnetic resonance close to lu = u^, which is governed by \/X. The expression for 
A, Eq. (IHUj) . suggests that the splitting can be tuned by the external magnetic field. The 
line width of those two resonances is defined by a' and /3, and the shape is a combination 
of the symmetric and antisymmetric Lorentzians as in Fig. HI Let us make estimates 
for a system shown in Fig. |6l The dimensions of the metallic links are here chosen as 
(0.5 X 0.05 X 0.01) /im with a Py load of the size (0.1 x 0.2 x 0.02) /im, for which we can 
estimate a resonance frequency in the range of 0.5GHz (/x ~ lOOGPa)^ and a coupling 
parameter g ~ 0.001 (free-standing metallic structures of such dimensions have already 
been realized^i). For intermetallic interfaces, the"effective field effect" due to 772 is very 
small, thus the phase shift of the magnetization is absent (Ty = 0). The amplitudes, given 
by Eqs. fl22|l23p . are large in the proximity of the ferromagnetic resonance (FMR). The 
magnetovibrational coupling splits the FMR peak by ^yguJ {Hq < N^Mg), which allows an 
electric detection of the mechanical motion excited by the spin transfer. Normalized voltages 
are plotted in Figs. [7] and |H] for purely ferromagnetic (mechanical motion is suppressed) 
compared to magneto-vibrational (MVR) resonances. 

Finally, we write the response function for arbitrary magnetization directions in Fig. (jS]): 



2eVm uj"^ — Wm + 2ia'uu^ + AF(ci;) ' 

where A is defined in Eq. fHOj) . As one can see, the form of the response function does not 
change. Eq. (1381) can be substituted into the first parts of Eqs. (I22|23l) in order to calculate 
the nonlinear response to an AC current. Similarly to Eq. flTB|) . the parameter Ty/ ^ 
Vi^x'y' + V^Nx' governs the balance between the symmetric and antisymmetric Lorentzians 




^7 sin 6 



(38) 



2eVm 



—irjiu + 'yMgTyf 
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Figure 7: Dependence of the DC component of the voltage on the frequency or the AC current bias 
for purely FMR (dashed line) and MVR (solid line) (w^ = Wg, a' = 0.02, /3/a; = 0.002, g = 0.001). 

composing each resonance (see Fig. H]). The magnetovibrational coupUng is described by A: 

Re A = ^ [w^ - {H,> + cot eiHy>>)u^/{N^,>Ms)] (40) 
Im A = gu'j [Hz'Nx"y"/Nx" + cot 9i{Hx" + HyiiNx"y" /N^")] (41) 

and 

g = sin^ OiM^VmN^'' {1/Ck cot(fcL)) (42) 
^ sin' e,M^VmN,.L/C = sin' 6, N,.{L/a)\V J Vi){M^ / fi), 

where Vm and Vi are the volumes of the magnetic load and the normal metal spring, re- 
spectively. L and a are the largest and the smallest dimensions of the normal metal links, 
k = Uc/ c, 9i is the angle between the equilibrium direction of the magnetization M2 and the 
current flow, and the external field is Hq = (i/^'/, i/^", i/^'). The parameter A is calculated 
in a reference frame x", y", z' shown in Fig. (|3]). The a;"-axis is perpendicular to the current 
flow and the axis z' (in general, this reference frame is different from the one used in Sec. 
IIIA, since the a;"-axis is not necessarily along the direction of rrii x 1112). 



17 




Figure 8: Dependence of the second harmonic of the voltage on the frequency of the AC current bias 
for purely FMR (dashed line) and MVR (solid line) (w^ = Wg, a' = 0.02, /3/a; = 0.002, g = 0.001). 

C. Mechanical torques due to absorption of transversely polarized currents 

Let us now apply a DC current bias to the system in Fig. O According to the magne- 
toelastic equations discussed above, the angular momentum of the spin polarized current 
is completely transformed into mechanical torques when the crystal and shape anisotropics 
are strong enough to prevent magnetization motion relative to the lattice. The mechanical 
torque then equals the spin-transfer torque: 

h 

Tst = -7— /r^im X (mfixod x m) ^ -7^0^71, (43) 

The Gilbert damping does not appear explicitly in this formula; however, it determines the 
time scale 1 / (ajMs) in which the system reaches quasi-equilibrium. The upper boundary for 
the spin-transfer-induced mechanical torque for a thin film of the size (20 x 200 x 200) nm 
(without crystal anisotropy) is defined by the maximum effective magnetic field due to the 
form anisotropy, Hejf = AttOMs, where 6 ~ 7r/20 ~ Ho/Mg corresponds to some small 
stationary deflection of the magnetization out of the plane of the film in the presence of the 
stabilizing external magnetic field Hq (Fig. [6]). This field generates a torque of 'jAnOMs ■ 
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MWNT MWNT. 




^ F2 



Figure 9: A spin-transfer nanomotor in which a metalHc wire (rotor) with strong spin- flip scattering 
(such as Pt) is grown inside a MWNT that is connected to two ferromagnetic contacts. 



MsVm/'^ ~ 10 ^^Nm, where Vm is the volume of the ferromagnet and Ms = 10^ A / m is 
the saturation magnetization for Py. Such torques are well above the sensitivity of existing 
NEMS oscillators.-^ 



IV. SPIN-TRANSFER NANOMOTOR 



Finally, we address the question how the torque can be transformed into a potentially use- 
ful rotary motion. It has been suggested in the literature to use carbon nanotubes as bearings 
for metallic nanowires.-"^ In Fig. [9l we propose a design of a spin-transfer nanomotor based 
on multi-wall carbon nanotube (MWNT) connected to two ferromagnetic electrodes (the 
torque doubles when the second electrode is ferromagnetic, but one ferromagnet is sufficient 
in principle). A metallic nanowire with strong spin-flip scattering nanowire is encapsulated 
by the MWNT. Pt would be a good choice (/f^* ~ 20 nm at 4.2 K^^ and /f^* ~ 1 nm at room 
temperature^^). FeCo nanowires that have already been grown inside nanotubes,— presum- 
ably have spin-flip diffusion lengths not much different from FeNi (/^^ ~ 5nm),— which 
makes this material also very suitable. The metallic nanowire should preferably be longer 
than the spin-diffusion length in order to achieve a complete angular momentum transfer. 
After burning off the outer shells over the platinum nanowire (in addition, the MWNT may 
also be pulled out to open the Pt wire), we force the current to flow through the metallic 
wire that serves as a spin-sink and a rotor. It has recently been calculated that the con- 
ductance may rapidly oscillate due to quantum interference effects as we change the overlap 
between two nanotube shells,—"— but disorder strongly enhance the intershell conductance, 
consistent with experiments.— 

The situation with not too high polarization of the MWNT connection to ferromagnet^ 
(P ~ 0.01) has tendency to improve with recent experiments reporting TMR = 5% (P ~ 
0.2).—"— Half-metallic contacts to nanotubes might lead to much higher values.— Adopting 
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the maximum current through a sine; le MWNT measured to date^^ J = 1 mA, we arrive at 
an estimate for the mechanical torque generated in the rotor T ~ Tq ~ 10^^^ Nm where we 
optimistically assume P ~ 1. This exceeds by many orders of magnitude the torques that 
can be induced by circularly polarized light (T ~ lO^^^Nm).^- The advantage of a MWNT 
bearing is that the friction force can be very small. The bearing may still get stuck at certain 
preferred positions that minimize the interaction energy of the sleeve and shaft. However, 
each layer is likely to have a different chirality, so the potential barrier (static friction force) 
hindering rotation should be very small. An upper boundary for the static friction force per 
unit area has been measured^^ to be 6.6 x 10^^'^ N / nm^ (but the actual value is possibly 
much smaller). For a nanotube radius of 2nm this corresponds to a static friction torque 
per unit area of less than 10~^^ Nm/ nm^. For the same nanotube radius, the overlap length 
between the rotor and the outer shells can be made at least 10 nm without the risk that 
the rotor gets stuck by the static friction. When the barrier to rotation is overcome, the 
static friction is taken over by the dynamic friction. The latter is much smaller initially, but 
increases proportional to the angular velocity.— 

The proposed motor could be useful in the next generation of synthetic nanometer-scale 
electromechanical systems. With an attached metal plate the rotor can serve as a mir- 
ror, with relevance to high-density switching devices.- The motor can find applications for 
inducing and detecting motion in microfluidics systems and biological systems. 



V. CONCLUSIONS AND OUTLOOK 



We find that the electric current-induced creation and detection of mechanical torques 
in magnetic nanostructures is possible by rectification technique suggested recently in Ref. 



9 and 



10 



We develop a detailed theory of the magnetization dynamics in the presence 
of electric currents relevant for these experiments that agrees with the very recent study 
by Kupferschmidt et al..^^ Subsequently, we predict that an alternating current can drive 
magnetovibrational dynamics, which can be read out by the generated DC voltage. Finally, 
we come to the conclusion that electric current-induced mechanical torques can create rotary 
motion. We propose a novel spin-transfer driven nanomotor based on integrating metallic 
nanowires with carbon nanotubes. 

In the study of the rectification effect, we find that the DC voltage originates from 
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two mechanisms: the rectification of the apphed AC current and the spin-pumping by the 
preces^ing ferromagnet. The second mechanism is important only in very thin ferromagnetic 

layers when the Gilbert damping is strongly enhanced (in accord with Ref. 35). It should 
be noted that the shape of DC voltage as a function of frequency is a symmetric Lorentzian 
for the spin-pumping mechanism. This can be distinguished from the voltage induced in the 
presence of noncollinear magnetic anisotropics or the "effective" spin-transfer field, which 
causes asymmetric line shapes. 

We generalized these results to treat magnetovibrationally coupled systems. The 
strongest coupling is achieved when the lowest mechanical mode is at resonance with the 
FMR frequency. In that case, a magnetopolariton is formed, and the Lorentzian shape of 
the DC voltage splits by an amount that is governed by the sensitivity of the mechanical 
system to external torques and by the magnetic anisotropics (without which there would be 
no magnetovibrational coupling). The technique based on the resonant magnetovibrational 
coupling can also be used to detect vibrations that are created externally. It can therefore be 
an alternative to the magnetomotive technique employed in fast transducers of mechanical 
motion.— 

We conclude that the functionalities of the spin-transfer torque, already used in appli- 
cations such as magnetic memories, can be extended by taking into account the coupling 
with the mechanical degrees of freedom. The experimental realization is a challenge, since 
free-standing metallic small structures on micro and nanoscale need to be fabricated and 
manipulated. 
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Appendix A: Effect of mixing conductance on the torque and magnetoresistance 



In this Appendix, we generalize the relations derived in Ref. |1J for the torques and 
resistances in a general N|F1|N|F2|N multilayered structures without bulk layer spin-flip 
scattering to the presence of an imaginary part of the mixing conductances or"effective fleld". 
We use magnetoelectronic circuit theory as formulated by Eqs. ( fTT][T2i) . Our result for the 
spin-transfer torque acting on the flrst ferromagnet (the torque on the second ferromagnet 
can be obtained by permutation of indexes 1 and 2, and by replacing / with — /) can be 
summarized by: 



rT.1 
J- c 



ST 



Ih 

Ye 



[?7smi X (mi X 1112) + r7/(m2 x mi)] 



(Al) 



where we introduced the spin-transfer efficiencies: 
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with a = cos 9, 4i?i(2) = 1/G'i(2)t + 1/^1(2)4 - 2/G'i(2)r-, 4i?i(2)_ = 1/G'i(2)t " 1/^1(2)4,, 

4/(^1(2) = 1/Gi(2)t + 1/^1(2)4, 2_Ri(2)r- = 1/Gl(2)r, 2/?^ = l/Cir + l/G2r and 6*1(2)1 = 

G'i(2)j/G'i(2)r; whcrc Gi(2)t and Gi(^2)i ai'e conductances of the left (right) ferromagnet in- 
cluding the left (right) normal layer (when necessary, the middle layer conductance can also 
be included in these conductances), Gi(2)t4, = G'i(2)r + ^G^i(2)i is the mixing conductance of 
the left (right) ferromagnet. 
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Figure 10: Enhancement of the torque by a nonzero imaginary part of the spin-mixing conductance; 
dashed hne for symmetric (P = R^G = 0.7, G/Gr = 0.6, Gj = 0.6) and soHd Une for asymmetric 



junctions (Pi = Ri^Gi = 0.6, Gi/Gir = Gi/Gir = 0.6, Gu = -G2t = 0.4, Gi = 2G2, P2 = 
^2-^2 = 0.13). 

The angular magnetoresistance reads: 



When the imaginary part of the mixing conductance is small, the absolute value of the 
torque Eq. (lAip and the magnetoresistance Eq. (1A2I) contain only second order corrections 
in Gii and G2i. In general, we conclude that the torque is enhanced when the imaginary 
part is not zero, as it is shown in Fig. [TUl This enhancement makes it impossible to have 
noncoUinear points of zero torque reported earlier— >^ (see Fig. [TOl) which can influence the 




^ [i?^_i?2 + -R2--^l] Rr -\- R2rGi^ -\- RirG2i + Ri_RlrR2rG2i -\- R2-RlrR2rGij^ 



(A2) 
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stable precessional states. 



Appendix B: Spin pumping and DC voltage in asymmetric N|F|N structures 

We consider here an asymmetric N1|F|N2 layered structure connected to two reservoirs 
and excited by rf magnetic fields. We are interested in the DC voltage that builds up under 
the condition of zero charge current. The system can be described by the generalized LLG 
equation and magnetoelectronic circuit theory: 

-7M X Heff + — M X — - -i— m X [(1,1 + 1,2) X m] (Bl) 



dt ' Ms \ dt J 2e14 

= (Gt + ) - f^f^) + (Gt - GM'^ - ■ m (B2) 

1.1(2) = m [{G, - ) - + {G, + ) - f^T^] 

- (2m X ^^(2) _ ^j^^ X mGr - (2m x p^^^) _ ^-^^^ q,^ (gg) 

where G^ and G^ describe the spin-dependent conductance limited by interface and bulk 
scattering of the ferromagnet, G^^ = Gr + iGi is the interface mixing conductance of the 
ferromagnet, the vector m = M/Mg is the direction of the magnetization, /iq and /x, are the 
chemical potential and spin-accumulation in the normal metals, respectively. 

For metallic interfaces, Gi usually amounts to only a few percent of Gr- We calculate 
here the linear response to the external rf magnetic field when Gi is disregarded: 

Xxx^ = (MJK)^ = ' ' . B4 

Xyx w = [My h^)^ = — — ^— — . B5 

where a is the bulk Gilbert damping, age = 7^^ [Gi + G2) / {2eMsV) is the extra damp- 
ing due to spin-emission, 1/Gi = 1/gi + l/2Gr, I/G2 = 1/(72 + l/2Gr with gi and g2 
being conductances of the R1|N1 and N2|R2 interfaces, respectively, a' = {a + ase)iNrc + 
Ny)/ {2^NxNy — N^y) (uotc that, when important, bulk scattering can be easily included 
into the conductances gi, g2 and Gr)-—'^ 

By solving Eqs. flBlHBSP to the second order in a small rf magnetic field the generated 
DC voltage can be expressed via the susceptibilities Eqs. (IBSp and (lB4p as: 

eV = 2nwPGr (- - -) Im (x^^xU) hl/M^, (B6) 
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where we introduced the effective polarization of the device P = — 1/G^/{1/G^ + 

1/^4 + 4/(71 + 4/(72)- The asymmetry of the interfaces R1|N1 and N2|R2 is seen to be crucial 
for the generation of a DC voltage in this geometry. Note that the dependence of the DC 
voltage on the frequency has a Lorentzian line shape as follows from Eqs. ( lB4tlB6[) . When 
P > 0, the positive sign of the voltage corresponds to the voltage applied to the junction 
with higher conductance (71(2) • 

The first order correction in the mixing conductance Gi leads to a slight change of the 
height of the Lorentzian and to a shift of the resonance frequency: 

,,2_^.^, in\Gl + Gl) G 

Note that the magnetovibrational coupling does not change the form Eq. (lB6p . when the 
proper susceptibilities Eqs. are substituted. The magnetovibrational coupling can thus 
be observed as a splitting of the Lorentzian peak due to spin pumping in an asymmetric 
N1|F|N2 structure. Generation of a DC voltage by magnetization precession of a single 
ferromagnetic layer has been suggested in Ref. |37j. However, those authors concentrate on 



the generation of a DC voltage due to spin-flip scattering in the ferromagnet. 

Appendix C: Spin pumping and rectification of AC currents in F|N|F|N structures 

It follows from the previous appendix on N1|F|N2 structures that a spin-coherent 
F1|N|F2|N structures in which the normal metals are identical and an extra ferromagnetic 
layer Fl (with fixed magnetization) causes the asymmetry, should generate a DC voltage 
as well. In such a structure, AC currents instead of the rf magnetic field can generate spin 
transfer torques when the magnetizations are non-collinear. Since the analytical expressions 
for arbitrary angles between the magnetizations are complex, we concentrate here, without 
loss of generality, on a 90 degree configuration. 

We assume a constant AC current bias on the system, thus forcing the extra charge current 
due to spin pumping to vanish. Following the derivation in Appendix B, and disregarding 
the imaginary parts of mixing conductances for both interfaces, we arrive at the following 
expressions for the linear response functions: 

.X (M in r^i YiM,{Ny + N,ya)-iujaf] 

XxiM - (M,//). - 2eK,^2(i + «se«.e)_^2^(i + ^2)+,^A' (CI) 
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where the anisotropic Gilbert dampings due to spin emission are a 



^h^G^ + G2)/i2eMsVj, a- = 7^' (^n + G2) / (2eM,Kn), A = 7M,(iV,(a + af ) + 
Ny{a + «-) + - OiV,,), l/d = (l^t + l/m)/4 + 1/2^, I/G2 = I/92 + 1/2^, 

l/G-^; = l/2Gr + 1/25(1,, with 5(1^,^ = gi^ + igu (however gu is neglected) and gi^^^ being 
the mixing and normal conductances of the R1|F interface, respectively. 

The second order analysis provides us with an expression for the DC voltage. After 
combining it with Eq. fj22|) . we arrive at the full expression for the DC voltage in terms of 
the susceptibilities Eqs. (lClllC2p . consisting of separate contributions due to rectification 
and spin pumping: 

- W-f-^''^^" + if, - S) -CX..X,.)/„VA«. (C3) 

where an effective conductance is 1/gi = {3/2gir + {1/Agi^ + l/Agi^)R^^/Ri) /4 and an ef- 
fective polarization of the device is P = (1 /4G| - 1 /4G^)/(1 /4G^ + 1 /AG^ + 1 /2gir + 1 /g2). 
The dependence of the first term in Eq. ( 1C3P on the frequency is in general a combination of 
the Lorentzian with an antisymmetric Lorentzian, as it is discussed in the Section IIIA. The 

second term becomes important when a'^^ > a and — ^7^^ > 1. When the damping 

due to spin emission is small or the asymmetry of the tri-layer weak, the second term in Eq. 
f ICSP can be disregarded. Its dependence on the frequency has a Lorentzian shape, see Eqs. 
( ]ClllC2p . In general, the signs of the first and second term in Eq. ( 1C3P can be opposite, 
thus possible suppressing the symmetric-Lorentzian part of DC voltage. 

The first order corrections in a small mixing conductances, Gi and gu lead to a change 
of the height of the symmetric and antisymmetric parts of Lorentzian and also to a shift of 
the resonance frequency, as it arises from the expressions for the susceptibilities: 

^ (,A-(M/n -J^ iu{r,2-r,ial^)+lMj:, 

[My/i)^ 2eK, co^il + a%-al- + 2k) - <(1 + a^) + luA ' 

XA^) - [M.JI)^ - 2ev;, + afa- + 2k) - <(1 + a^) + .a;A' ^^^^ 
^2 \ G- 

where k = 1 y (;^~^ + ;^)^' = ^^^^v + ^^2/") + ^^^^^v " ^v^) ^"^^ = 
r]^{N,y + N,a) + ^2 C^x - N^yO). 
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Finally, we present our results in the presence of the magnetovibrational coupling {e.g. 
see Fig. |6]) keeping only the dominant terms in small Gilbert damping: 

^ sin 6*7 

iu{rj2 - Viatf + ^(1 - N,jN,)Fiu)) + jM,T, 

9->o h sin 6'~f ri2iuj + 'jMsT^ 
^ 2eVm a;2 - + iuA + AF{uj) ' 

h sin 6'j 

-^UJ{'n^ + rj,K + r/^Q + 7^,1, [1 + ^ {HjM.Nl) F{u)] 

^ uj\l + 2ti)-ujl + iujA + KF{uj) ^ ' 

g^o h sin 6'-f —rjiiu + '-fMsTy 

2eVm — Wm + ii^A. + AF{uj) 
These susceptibilities can be used in Eq. ( 1C3P in order to calculate the generated DC 
voltage. In the regime of resonant magnetovibrational coupling, the second (spin emission) 
term in Eq. (1030 is comparable with the first one only when a**^ ~ P/^m- The second term 
corresponds to two symmetric Lorentzian peaks split by Re A. 

Spin pumping therefore can play an important role when the enhanced interface damp- 
ing is comparable to the bulk Gilbert damping and, in the regime of magnetovibrational 
coupling, the mechanical damping. 
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